General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



nu/TMM 82-4 


SEHt-AMNUAL Pnoeffiss Report on 

THREE DIHEHSrOW. aOH FIEIS IHS1!£ OJHPKSaW 
ROTOR^ INCUOIIK BLASE BOtWAW LAYERS 


J. N. 6AUCS, n. POUAGARE^ AND B. LAI^INARAYMA 


C14SA-CB-169120) THEEE DlflEBSIOMAi ELOi 

field inside COHPBESSOB fiCXOB* INCLOD1I6 

BLADE BOUBDABT LAIEBS Se»ia»Bttal Progress 

Beoort iPeniisyleaiiia state Ociir.) 78 p 

BC A03/HF AOl CSCL 20D G3/3M 


M82-27686 


(jBclas 

25219 


NASA Grant NSG 5266 

National Aeronautics % Space Administration 
Lewis Research Center 


i 

i 

Turbomachinery Lak>ratory 
Department of Aerospace Engineer ins 
The Pennsylvania State University 
University Parr^ PA 16^ 


PSU>Turbo>82-4 


S«ii-Aiumal Progress Report 
on 

THREE DDfENSIOHAL FLOW FIELD INSIDE C(»mESS(» ROTOR 
PASSAGES, INCLUDING BLADE BOUNDARY LAYERS 


J. M. Galmes, M. Pouagsre, ai^ B. Lalcshalnarayaiia 


to 


NASA Lewis Research Crater 
Project Monitor: Dr. P. Sockol 


Turbonachlnery Laboratory 
Department of Aerospace Engineering 
The Pennsylvania State University 
University Park, PA 16802 


June 1982 



11 


PREFACE 


The progress of research on "Three Dlaenslonal Flow Field Inside a 
Oppressor Rotor Blade Passage, Including Blade Boundary Layers" (KASA Grant 
NSG 3266) for the slx-9K>nth period ending June 30, 19B2, Is briefly reported 
here. Two papers were presented and published during this period. These 
are listed in section S of the report. 


B. Laksholnarayana 
Principal Investigator 
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1. TURBULENCE MODELLING: REYNOLDS STRESS MCH)EL (RSM) 


A literature survey on both the analytical and the experimental work 
on effects of curvature and/or rotation was given In ref. 1. Some of the 
Important conclusions of that survey were: 

1. Only few calculations are available for the prediction of the 
three dimensional boundary layer In rotating frames. 

2. No complete Reynolds stress model Is available for rotating turbulent 
flows. Very few attenq>ts have been made to account for the rotation effects 
In the k-£ model and those are not based on a logical analysis. 

3. There are very few detailed measurements providing Information on 
the effects of both the Rossby and the Richardson number on turbulence. 

It was emphasized In ref. 1 that new experimental results would be of 
great Interest, particularly if the rotation effect can be Isolated from the 
other effects. It was also noted tha t a major effort should be given to 
the analysis of the dissipation rate equation and to the Reynolds stress 
equations . 

In the previous report tll» modelling of the rotation effect and the 
low Reynolds number effect In the dissipation and the kinetic energy equa- 
tions were described, with major emphasis on the equation for the dissipation 
rate. It was also noticed that. In the case of a k-c model, the assumption 
of the existence of an eddy viscosity concept could not account for the 
anisotropy of the turbulence which exists In the boundary layer around a 
blade. The best way to avoid this problem would be to solve the complete 
set of Reynolds stresses equations along with the momentum equations. 

However, the calculation procedures for three-dimensional viscous flows 
require large memory storage and large computer CPU t.'.me to solve the three 
momentum and continuity equations. The resolution of tre six Reynolds stress 
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equations sintultaneously with the three monentua equations, the continuity 
equation, end an equation for the dissipation rate Is a very cosq>lex problea 
and can still be considered as an unresolved one. It Is hoped that the 
rapid progress In computers and also the progress made In developing calcula- 
tion algorithms for three-dimensional viscous flows will allow such calcula- 
tions to be done In the very near future. 

A two-equation model Is a compromise between a full Reynolds stress 
model and the empirical models which fall to rep «.v sent the turbulence 
properties in such complex flows as rotor blade boundary layer. However, 
the anisotropy of turbulence which exists In flows cannot be represented 
by the usual Isotropic eddy viscosity formulation. In Its more general 
form, the eddy viscosity Is a fourth order tensor which depends on many 
parameters such as the Reynolds tensor, the strain tensor, the curvature 
and the rotation. It Is therefore, a very complex formulation, which evidently 
represents very complex phenomenon. At this stage of the discussion It 
appears that the eddy viscosity tensor could be as difficult to handle as 
a Reynolds stress model. For the cases which are of Interest to us, e.g. 
blade boundary layers, an algebraic model of the Reynolds stresses based on 
the complete model, together with the dissipation equation might be a 
compromise between a detailed description of the turbulent stresses (full 
Reynolds stress model) and the very crude eddy viscosity hypothesis. 

It Is not necessary to solve the complete set of equations for the 
Reynolds stresses, it Is, however, necessary to analyze the Reynolds stress 
equation to find out fetich terms have to be modelled or neglected and then 
derive some propositions to model the remaining terms. 
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1 . 1 Reynolds Stress Equation 

The equations of the mean and turbulent quantities for an Incompressible 
flow In the generalized tensor formulation are given In ref. 1. The Reynolds 
stress equation which Is of Interest In the discussion that follows Is 
given by: 


uJ),j - sJ + p uJu'u’J - ujrJ - 


+ - Pu;u'^(n ^ +2n^e^ J - pu>'J(u, , + ^ 

i,J k k,j 1 1 k,j kpj k 1,J Ipj 


u ' - u ' 

k l.J "i kj 


Ik 


2y S 


Ik 


( 1 ) 


It Is evident that the Reynolds stresses appear to be affected explicitly by 
the Coriolis forces, but they are also Implicitly affected by the rotation 
through the triple velocity correlations, the pressure velocity correlations, 
the pressure strain correlation, the production by the stresses themselves, 
and the dissipation. 

In the case of a Reynolds stress equation model, the second order terms 
are handled exactly, then only dissipation, pressure strain correlation and 
diffusion terms need to be analyzed and modelled. 

These quantities must be represented as empirical functions of the 
mean velocities, Reynolds stresses, and their derivatives and the rotation 
and curvature. These terms will be analyzed and modelled separately. 


1.2 Pressure-Strain Correlation 




i.j k 


6 ^ + p ' u* 


kJ,J ^-P^,k^P\,l 


Following Chou [2], the explicit appearance of the pressure may be eliminated 
by taking the divergence of the equation for fluctuating velocity u^, thus 
obtaining a Poisson equation for the fluctuation p'. 
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Th« fluctuating Navier-Stokea equation la: 


(puJ)+(pU|^u’’^+ pujU*^ + 


let ua then take the divergence of equation 3 


(3) 


g‘“(pu’^).^ + »*”(P0(U’' + puju' + pu'u*' 


Puiu'F).„ + 


2g*”'(et^^pn'u’'),_j • 


and with the continuity equation: ; uJ^-O 

we nay re-wrlte thla equation. 

a) the time and apace derivative are consButableg then the time derivative 
dlaappeara 

b) the vlacoua tern alao dlaappeara. 


which la the Polaaon equation for the preaaure fluctuation. 

Then following Chou, the preaaure fluctuation may be expreaaed In 
the following form for a poaltion x In the flow. 

dV 

^ ^ ///™i ^ 

dV' 

+ :f /// , (5? + 

2TT •''-’vol ‘ l,r, pr 1 l.n. i |?i 


1 
3p; 




*1 1 If I 
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(4) 


(5) 


where tema with and without an Indice 1 relate to valuea at and x, 
reapectlvely (the integration being carried out over apace) (C - Xj^ - x) . 
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An cxprAtslon for the pressure strsln correletlcm which eppeers In the 
equstlon of Reynolds stresses cen be derived es follows. We celculete the 
correletion et point x, then ^ end ^ ere Independent of point end 
mey be cwsldered es constents for the volume end surfece Integrels. 

Th«i we mey write: 




^ * ♦tk.2 * 





ES 


♦ki.l * *kl.2 * Set ■ »kl'*!.l + + 'J’ 


( 6 ) 


with: 


/ - A ff r (u'\trv ^JlL 

^ k,l 4 tt -'JJvol '‘"l “l ^’nr^ ’k |j|i 


//J , {u® + e" ul^ t, 

k,2 2 it ^^^vol ^ l.rj^ pr 1’ l.m^^ t,k 


dVj^ 




^k Axp ^jjljj 3u^ “l.k *^1 l,k an^^ 


(7) 


Equetlon 6 suggests thet there ere three distinct kinds of Interectlon giving 
rise to the pressure strein correletion; one involving fluctuetlng quentltles 

I 

. , snother erlslng from the presence of "extemel effects" such es meen 

Kf 1 

I 

strein rete end rotetlon <1* 2’ * surfece integrel which 

will be negligible ewey from the vicinity of e solid boundery [3]. 

Some of the ptoposels for closing the Reynolds stress equetlon heve 
essumed thet ^ Is the only slgnlflcent contribution to pu^ ^ [4-6]. 

However, Reynolds [7] hes shown thet prediction of e renge of even homogeneous 
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frM turbulent flout deiaendt the Inclueton of sMen etrein rates In the 
preeture-strain terms. Moreover, both Townsend (8l and Crow [9] have shown 
that under conditions of rapid distortion the effect of outweighs 

th« of 

It is then necessary to model 2 *• i* 

1.2.1 Modelling of Term" 

Our practice in siimilatlng 2 2 direction 

from Rotta's analysis and more recently Launder, Reece, Rodl (3] and Lumley 
[11]. It is easy to show that the rapid term can be written as (Appendix A) 


♦lk,2 * *kl,2 


< <*Lk * Vi> 


where 


•*m _ rJ 


+ e“ 
’r pr 


•mlk " " 2 it ^“1 “1 ^’mj^k 


J!!L 

Hill 


( 8 ) 


Equation 8 is a rigorous consequence of equation 6 when all second deriva- 
tives of the mean velocity are negligible and the turbulence field is 
homogeneous. It is of course only approximately true in more general flows. 

In the case of a cartesian coordinate system, Rotta [10] has commented 
that the fourth order tensor should satisfy the following conditions: 


•>™-*try= ‘ ‘mlk ' *tkl. 

inccmpressibllity: ® 

normalization ' 


rikk r i 


(9a) 

(9b) 

(9c) 


rmlk - «tl - «rk «,il>T5 


isotropic turbulence: a 


(9d) 
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Th« spectrum of the Rsynolds strses cooplstsly dotendnst th« form of a .. . 

tmiK 

If va aaauma chat tha apactrum haa an aqulllbrium form, wa may axpact tha 
form of chla fourth ordar canaor to ha axpraaalbla In carma of tha Raynolda 
acraaaaa. Thara ara 10 llnaarly Indapandanc coi^inaclOTia of tha anlaocropy 
canaor which satlafy tha aymaMtry emdittona, incompraaalbllity, normalisa- 
tion, Mid laotropy condition. In practlca, tha flrat two llnaar tarma ara 
uaad [11]. Howavar, Lumlay haa shown chat this approximation doss not 
satisfy tha raalixablllty condition. It la nacaaaary to go at Isaac to 
sacond ordar tarma bafora this can ba achlavad. A nodal la said to ba 
raalltabla If It guarantaea that quant It las which should ba non-nagatlva 
(Ilka varlancas) will ramaln non-nagatlva, and Chat corralatlon coafflciants 
will navar axcaad unity In absoluta valua. 


Lumlay* a Hodal fill 


rmlk rl ak rl w rl ok rm Ik rk ml 

+ c(b 5 . + b .6, b,. - .b. ) (10) 

m Ik im rk rk la Ik cm 3 rl ak 3 rl ak 


with 


u'u’ , 

b j ■ ~^T“ " T 5 4 and a .. 
rl 2 3 rl rmlk 
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2kl 


rmlk 

and only ona constant to adjust. 


Laundar, Raaca and ltodl*a Modal [3] 

Thesa muchors approxlaaCad the canaor by a llnaar combination of 

Ryanolda atrassas. Thalr mathod la vary similar to Lumlay* a; tha synMCry 
constraints l^ily that tha fourth ordar tanaor may ba vrlttan aa: 
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a ■ E6, u'u' + F(6, . u'u’ + 6 . u’u,' + 6 . u'u' + 6 u’uj) + C_5 . u,’u' 
rmlk kn r 1 ki r m ml r k kr 1 m mr i k 2 ri k m 


+ G 5 It + SC'S, .5 +5 ^5, )k 

rl km kl mr ml kr 


( 11 ) 


frtiere E, F, G, are constants. The application of equations 9b and 9c 
enables four of these constants to be expressed In terms of the fifth: in 

terms of 
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2OC2 + 6 
55 


2 


In cartesian coordinates, equation 8 may be then written in the following 
compact form: (Launder, Reece, Rodl [11]; Coustelx, Aupolx [13]) 


'*’ik,2 \i,2 


C, + a - 30C- - 2 3U* 3U* 

_2 (p* _ 2 ^ , 2 __ki 

11 i^ik 3 ik*^ ^ 55 

^ K - f 


( 12 ) 
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i"j ^ 3^ 


3U* 3U. 

+ e jjP 
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It is interesting to note that equation 12 has been derived mainly from 
kinematic arguments, for the case of isotropic turbulence subjected to 
sudden distortion, equation 12 reduces to: 

Ic^ 

+ " 5 — ) Irrespective to the constant C_ 

oX^ Z 

which is the exact result derived by Crow [9], 


*kl>2 


3U. 
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However, It appears that equation 11 together with relations 9 cannot 
be generalized to any kind of coordinate system. In fact, if we go back to 

u* 

equation 8 it is easy to show that the fourth order a^^^^ does not satisfy 
conditions 9b and 9c. A detailed wialysls is then necessary to derive 
a general formulation for relations. While we await the final resolution of 
this problem (the study is under development), it is Interesting to point 
out some features of this problem. 

It appears from equations 10 and 12 that the "rapid term" of the 
pressure strain correlation might be expressed in tems of the Reynolds 
stresses or the anisotropy tensor and what %fe will call here the "external 
effects", e.g. mean strain rate, rotation. Moreover, one of the essential 
properties of the pressure strain correlation is its character of redistrlbu- 
tlvity (e.g. this term vanishes under contraction of indices). It also 
turns out, looking at equation 12, that the "rapid part" of the pressure 
strain correlation, is a turbulence production like term, e.g. the term 
may be approximated by a similar expression to the production by the 
turbulence through the "external effects." 

Then a simple way to express the "rapid term" may be given as follows: 


♦lk.2 * \l.2 ■ * 5 


(13) 


where 


P* 


n*. ■ K.i * "*.j 

and 


In order to be sure that the expression will behave properly in all coordinate 
systems, it is necessary that the constant be a function of the invariants 
of the anisotropy tensor, and the local Reynolds number of the turbulence. 
These Invariants are defined as: 
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II • 



III 




- ja 
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(W 


A formulation of the coefficient C is under survey, some directions towards 

P 

this end have been proposed by Lumley [11]. 

According to equation 13, 2 2 isotroplee the turbu- 

lence production tensor. However, it is interesting to note that if the 
production due to the mean strain is following this trend, it is not true 
for the terms involving the Coriolis force effects. 


1.2.2 Modelling of "Return to Isotropy” 

Let us first present some basic features about this term: 


4 > 


<P 


ik,l 

kl.l 


— fff (u' 

4ir 1 1 ’im l,k 


dV 


Atr ‘^•l^-f^vol ^“1 “1 ^*Jm“k,i 


dV* 


(15) 


This term involves only the fluctuating quantities and Is responsible for 
the return of anisotropic turbulence to Isotropy, in the absence of other 
disturbing effects. Rotta (10] pointed out that is a symmetric 

tensor with zero trace, which vanishes if the turbulence Is Isotropic. Since 
this term acts to Interchange energy among the conq>onents when the turbulence 
is anisotropic, and vanishes when it is isotropic, it Is natural to express 
it in terms of the anisotropy tensor of the turbulence, which Is also a 
symmetric second rank tensor which vanishes if the turbulence is isotropic. 
It is also easy to show that, for a flow without mean strain and rotation 
and initially non isotropic the Reynolds stress equation may be written 
as (cartesian tensors) 
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2 

— — •(fc +<t> — — 6 e 

3t ^ik,l \i,l 3 ivr 

which leads to 




3t 


" ^ik.l \l.l 


Therefore, as the flow must return to an Isotropic state, the following 
assumption Is made: 


'^lk,l \l,l 




(16) 


where — Is the rate of time at which the energy Is mainly transmitted from 
big eddies to dissipative ones. The value of the constant has been originally 
set to about 1.4 by Rotta. Different authors have then used a value of 
of about 1.5. 

According to equation 16, the sign of ^ ^ Is always such as 

to promote a change towards Isotropy, Its magnitude being proportional to 
the local level of anisotropy. This term Is called "return to Isotropy" 
term, and this model gives acceptable behavior in most situations of practical 
Importance. However, It does not take Into accotmt the variation of the co- 
efficient with the local Reynolds nuviber of turbulence, or the variation 

with anisotropy. It is found experimentally that is larger when the 
turbulence Is more anisotropic. Consequently, the coefficient appears 
to be a function of the Reynolds number and of the anisotropy. This 
anisotropy Is generally created by "external effects", such as shape of the 
boundaries, rotation, mean strain. In fact, these effects mainly Influence 
the time of return to Isotropy. Then we may write that the tensor 

is a functional of the anisotropic tensor, but also depends 
on the Reynolds nun4>er, the "external effects" explicitly. In fact, we 
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may divide these effects in three main categories: 

the gross turbulent structure effects: through anisotropic tensor 

the dissipative effects: through Reynolds number of turbulence 

the "external effects": through mean strain, rotation, curvatures 

It is shown in Lumley [llj that, even though in many ^lows the turbulence 
does not become sufficiently anisotropic to make realizability a necessary 
'ondltlon, there are, however, some situations where the *'.urbulence does 
become nearly two dimensional. The constant should then reduce to 1. 

It is clear that in order to be sure that the expression will behave properly 
in all corodlnate systems and situations, it is necessary that be a function 
of the invariants of the anisotropy tensor, the Reynolds number and the 
"external effects". 

Lumley [11] proposed a formulation of the constant C^. In his paper 
the constant is 2C^. Then the formulation is: 


2C 


^ - 2 + F(Rp, II, III) (1/9 + 3III + II) 


(17) 


where Rj^ ■ q /9ev and II and III are the Invariant of the anisotropic 
tensor. 

F is then expressed as 

F - exp(-7.77/Rj^'){72/Rj'^^ + 80.1 Ln(l + 62.4 (-II + 2.3 III)}} (18) 


Equation 18 is simply an interpolation formula to connect the known value 
In a certain number of experiment. However, this proposition gives certain 
directions which could be followed to establish a formulation of the constant 

''i- 

Cj . II, III, (19) 

where Is the local Reynolds number and R^^ is the local Richardson 
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nundier of the rotation. The fonnilatlon of the quantity la under survey. 
It la not known now whether all these effects are really necessary to be 
taken In account. However, It Is certain that the Reynolds nund>er Is to 
be Important In the evolution of the coefficient C^. At this stage arguments 
concerning the expansion of In a power series remains one of the 
directions to follow to derive an expression for the coefficient. 


1.2.3 Modelling Near the Wall 

Near the wall, the surface Integral In equation 5 appears to 
become Important in ma^ltude and must not be neglected. In fact, Bradshaw [14] 
demonstrated that the surface Integrals and would make a significant 
contribution to the pressure-strain correlation as long as the typical size 
of the energy containing eddies was of the same order as the distance from 
the wall. As he remarked, this condition is always satisfied In near wall 
flows. 

If we refer to Launder, Reece, Rod! [3], for a plane wall (with x^ 
normal to the surface) the pressure strain correlation may be recast In the 
following form, from which the surface Integral Is eliminated: 




vol 


{ 




3u; 

^3Uj 


1 ^ 1 r 1 1 


1 

T i 



[|x - y| lx - y*lJ 


dv 

( 20 ) 


where y* is the Image of the point Y in regard to th>^ wall. 

This form again suggests that there should be two contributions to the 
near-wall effect corresponding to the reflecting wall Influence of (^ . and 

XK f X 

<^lk 2 * Therefore, they proposed a formulation of the wall proximity effect 
on the pressure strain term as follows: 
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(♦lit + ■ '<=1 ! - 1 * li- ('■S * C f «“ 

n 4L 

where I denotes the length scale of the energy containing eddies, thus 

3/2 


fB - 


CXo 


ds 


In the general case, the two surface Integrals are: 

i- - f u' - I 

Hk 4 ttp 3nj^ l,k *^1 l,k 3n^ ' 


5 f {-^ ^ u’ - 77^ da 

^kl 4TTpJsi|[^|| 3n^ k,l ^1 k,l 3n^ ' 1 


(22a) 


(22b) 


where the quantities with subscript 1 correspond to point p^ and others 
to point p. 3/3n^ Is the normal derivative and dS^ Is a surface element. 

The value of this Integral Is to be evaluated at point p. Following Raj [15], 
a series expansion would Indicate that, to a first order approximation, 
equation 22a can be written as; 

4iTl|f|I 

It Is clear that near the wall the devlatorlc of the Reynolds stress and 

the "external effects" should have some contribution to the surface Integrals. 

Unfortunately, the experimental results very near the wall are extremely 

rare, particularly on curved surfaces, and It is therefore very difficult 

to derive a specific modelling of the surface integral In the case of 

highly curved flows. However, for flows over mildly curved surfaces we 

may hope that the only important length scales are the local turbulen t 
k3/2 

length scale — ~ — and the normal distance to the wall y. Moreover, it Is 
also hoped that the predictions will be less sensitive to assumptions made 
to third order correlations, handling the second order ones exactly. Then, 
following an analysis close to Launder, Reece and Rodl*s one, the effect of 


f 

l.k ■ 


jdeviatoric, external effects} (23) 
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the surface Integral is modelled through the Introduction of a function 
of the turbulence length scale and the normal distance to the wall y, in 
the general model for (p\i* . + p\i* .) . 


P’u! . +Vu^ 




where F( r) must -reduce to one away from the wall. 


1.2.4 Summary 

A model for the pressure strain correlation has been proposed, 

which Includes the effects of rotation and curvature explicitly. However, 

£ 

this model is still under study as three functions C, , C and F( 7 ) have 

1 P nj^x* 

to be defined. Some directions have been given, but further w>rk is needed. 

The pressure strain correlation model is presented above Equation 24). 

The coefficients C, and C are functions of the local Reynolds nuidier, the 
1 P 

Richardson number and the invariants of the devlatoric of the Re 3 molds stress 

I 

tensor, while F( r) is functlcxi of the length scale and the normal distance 

n^x* 

to the wall. 


1.3 Dissipative Teras and Diffusion 

Referring to equation 1 for Reynolds stress, two kinds of terms remain 
to be modelled — the diffusion term, which included viscous terms as well 
as pressure velocity correlations and triple velocity correlation, and 
the dissipation terms. Let us examine the dissipative terms. 
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Dissipative T^niw 

D - -fJu* - rju* where F . - U(u . + u )- 2vS . (25) 

k l,j 1 k,J IJ 1,J j,l IJ 

In Incompressible flows, only Che fluctuating quantities remain In equation 25 
then. 


Terms like a/'^u' . and s’'^u’ have been derived in the previous report and 

IC 1 i^ij 


it has been shown that we may write: 


®k i.j 3v ‘ Ik 1 e ik^ 


(27a) 


l*lk**2 '^I’lkl 


(27b) 


Therefore equation 26 may be approximated by 


® f '^k + • ■'i ‘ikl 


(28) 


A dimensional analysis, using similar approximations as the ones used In 

the previous report shows that the term Involving the mean strain Is of 
- 1/2 

order , and therefore Is negligible In high Reynolds nxuber flows. In 

that case the decay rate D reduces to: 

D - - f S,, pc (29) 


which is the form proposed by Launder, Reece and Rodl [3] who assume the 
dissipative motions to be isotropic. Several experimental studies have 
shown that turbulence does not remain locally isotropic In the presence of 
strong strain fields [8, 16]. Nevertheless, equation 29 seems to be enough 
for most of the flows studied by Launder et al. In the case of low Reynolds 
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nua^«r flows, £h« snsrgy containing and diaaipatlon ranga of motions ovarlap 
and.cha dissipation rats Is than connonly approxlmatad as [10] 


“i"k 


(30) 


Equations 29 and 30 have lad a nuobsr of workars to proposa that In ganaral 
the correlation may be approximated as: 


- - 2 pe {(1 - 27^ ^s^ 


(31) 


where is a function of the turbulent Reynolds number ^ * 

A detailed analysis of equation 31 has been proposed by Hanjallc and 
Launder [17]. However, equation 30 is not asymptotically valid. In fact, 
as the Reynolds nundier Increases the equatlm 30 should reach continuously 
the form obtained for very high Reynolds number (equation 29) . It seems that 
equation 28 should not present this problem, therefore we may re-write 
equation 28, Including the Reynolds nunber of the turbulence 




(32) 


where a is a constant of order 1. Some numerical calculations are needed 
to derive an "exact** value for the constant a. 

Diffusion Terms 

There are two kinds of "diffusion" terms, the viscous ones and 
the third order correlation ones. The diffusion terms In equation 1 are: 


triple correlation 


viscous 



18 


A dlacnsional analyili, Again shows that Che viscous tsxas are of order 

R ^ and smaller than other terms. Therefore, it is evident that these ten» 
e 

will be important only for low Reynolds nunR>er flows, very near a wall. The 
triple correlation terms need to be modelled so that the complete set of 
Reynolds stress equations is closed. This work is presently under survey. 


2. TORBULEKCE MODELLING: k-e MODEL 


2.1 Introduction 

An annlyalt of the kinetic energy end the dlselpetlon equations have 
been presented In the previous report [1]. At thet tlae, we sude the following 
remarks . 

The rotation does not appear explicitly In Che kinetic energy equation, 
but is present In the dissipation equatlcm. However, In high Reynolds 
nuad>er flows, we have shown that this explicit effect was negligible. At 
low Reynolds nu8d>er (the Reynol<5s number Is based on the following turbulence 
characteristics: k, e and Is v' ht • R^) , the dissipation becomes non- 

isotroplc snd both Che effects of Reynolds ntssber and Richardson nundter nay 
be important In this case. But Che main remark Is that the effects of rota- 
tion are more iiq)ortsnC on the production of turbulence Implied by the 
interaction of the Reynolds stresses and Che "external effects", than on 
Che others. In fact. If we refer to the previous chapter It Is easy to see 
that Che Reynolds stresses nay be greatly affected by Che rotation while 
the kinetic energy and Its dissipation rate are not so much affected. Then, 

It appears that in any cwo-eqtiatlon turbulence model combined with an eddy 
viscosity hypothesis, the first priority Is to control Che calculation of 
Che eddy viscosity coefficient. In fact. In most of the %K>rk done, the 
main hypothesis Co derive Che viscosity law assumes Chat Che turbulent viscosity 
is isotropic. For slisple shear flows, this gives adequate results at very 
low cost. Ho%rever, for three dimensional flows and particularly for some 
boundary layer flows the velocity vector u and 7u are not aligned in 
general. Therefore the isotropic eddy viscosity Is not adapted tc predict 
Che behavior of the Reynolds stresses. Solving the complete set of Reynolds 
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•treat cquatlone would avoid thia difficult problas, but it la a traaandoua 
work which involvaa vary inportant nuaarical achtMa* which arc only in thair 
•arly atagaa of davalopiaant. On tha other hand, the eddy viacoaity ia 
a fourth order tanaor in general, however, in laany real flow aituationa, 
there »ly two directlona in the flowa which are of equal inportance 
(it ia the caae for blade boundary layera). Therefore, we are able to 
derive different eddy viacoaitiea to take account for the aniaotropy of the 
flow. More generaiJy, it would be intereating to cleaaify the flowa that 
the engineer ia encountering in which could uae auch an eddy viacoaity 
model. Thia cvi be related to "Zonal Modelling" referred to by Kline during 
the laat Stanford Conference on coaqilex turbulent flowa. 

To be able to check oobm of the aaauoptlona made prevloualy, particularly 
for low Reynolda nuaber modelling, a coaq>uter code waa developed baaed on 
the Patenkar-Spalding procedure [18], which aolvea the parabolic two- 
dinenaional and quaai-three dimanaional tranaport equation for the velocity, 
the kinetic energy ctd the diaaipatlon. Chapter 2 ia divided into two parta. 
In the firat part, we preaent aoma aimple approach to the modelling ^he 

different atreaaea. In the aecond part, aome calculationa of a boundary 
layer developing on a rotating and non- rotating cylinder ia preaented. The 
non-rotating cylinder calculation haa been done to ahow the effecta of the 
correctlona due to Seynolda number end Richardaon nuaiber of rotation in the 
dlealpation equation. A brief erplanation of the Patankar-Spalding ntuaerlcal 
method ia given in Appendix B. 
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2 .2 AlMbf ic Mod«llin£ of the R«ynold« Stw»€i 

Th« Reynolds stress equation presented In Chapter 1 (equation 1} , 
toRSther with the essui^tions made, aey be written as: 

(» jijs’ u^).j ^ +«3pi) 

(33) 

where the left h«id side is the convective like tarn, the first group of the 
right hand side is the diffusion like tern «id C(uj^u^) is the source-sink ter«. 



The kinetic energy equation may be deduced from equation 33. 

(pk) + (okU''>.j - -(p^ * Oku’^ - ujF^^) ^(k) (35) 

1 Ik --r— r 

where C(k) is the source-sink term and k ■ g 


- P - e (36) 

P 

The transport end diffusion terms are treated by the technique known as 
;ilgebraic stress modelling (19). The net transport of u|u^ is assumed to 
be locally proportional to the net transport of k, the coefficient of 
proportionality taken as u’u^/k. Thus for incompressible flows. 
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(u|u^)+(u^u^ ~ • ~\~ '*' " Dif(lc)] (37) 

where Dlf denotes diffusion terms. 

This trestnent appears to be reasonably accurate for thin shear layers 
except near an axis of synaaetry. However, the turbulent boundary layers In 
a rotor are far to be axlsynoetrlc and more over If any syn&etry exists It 
must be a synnetry In regard to the axis of rotation which Is never in 
"contact" with the fluid. Therefore, substitution of the model assumptions 
yields the following algebraic equation for the Reynolds stresses: 



(38) 


r 4 ^ 

The form of the wall-damping function F[ rj miist now be specified, as well 

n^^x^ 

as the values of the coefficients C and C, , and a. This Is presently under 

P 1 

survey, however, we still may discuss equation 38. This equation is basically 
a non-linear algebraic system of six equations. In fact, which have to be 
solved simultaneously. The six equations plus a formulation for the dissipa- 
tion rate, will form a closed system, providing we know the mean flow. 

It is Interesting to note that several authors have treated similar 
equations as a linear system, provided that the kinetic energy k snd dissipa- 
tion rate c are calculated with two differential equations, and most of all 
the computations are restricted to local equilibrium turbulence for which 
P/e is unity. 

In this case, equation 38 provides mainly some indications on the 
anisotropy of the flow and then the quantitative values are given through the 
calculation of k and £. It is therefore not necessary to define any eddy 
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viscosity formulation. If we want to solve this system still linearly, when 

P/c is a function of the space. It is necessary to uke an assumption on 

the denominator in equation 38; then provided that the evolution of the 

stresses in the main direction of the flow is not very important, which is 

a good approximation in bomdary layers flows, we may calculate the stresses, 

supposing that P + e[c.F ( — - l] is taken at the previous step, in a 

' 1 J 

parabolic marching method, where all the quantities are known. 

On the other hand, it is also possible to solve the system of six non-- 
linear equations plus the equation for the dissipation simultaneously. The 
number of equations to solve is then reduced to seven. However, it is 
well known that non-linear systems are very much sensitive to the initial 
conditions we must give in order to iterate and converge toward the solution. 
But still, in the case of parabolic type flows, the evolutions in the stream- 
wise direction are very much less important than in the normal plane, so 
that it is possible to use the values at the previous step as a first guess, 
and as the changes should be small, the process should converge rapidly. 


2.3 Calculation of the Boundary Layer Over an Axial Cylinder 

The governing equations in the axial and circumferential directions for 
a uniform property, axisyanetric turbulent boundary layer flow may be written: 


M ^ 




9x r 3y 


- pu'v' 






1 

r 



3W/r 

3y 



which together with the continuity equation 


3rU 5rV 
9x 3y 


0 


(39) 



24 


and of radial equilibrium of the mean motion 



Constitute a closed set if we define a model for the two Reynolds 
stresses puv and pw'v’ . The independent variables x and y are respectively 
the axial and radial direction. The corresponding velocities are U and V. 

W denotes the circumferential velocity and is zero for non rotating cylinder. 
All symbols are defined In the nomenclature. The effective viscosity of 
this flow may be taken as the sum of the laminar and turbulent contributions, 
i.e. 


^eff " ^ ^ 


Provided that the two Reynolds stresses are respectively aligned with the 
correspondent mean strain. The turbulent viscosity is obtained from the 
solution of the following diffeirential equations. 

Turbulent kinetic energy k: 


3r y + — 


r. 3k . 3k 1 ' 

37 


3y 


+ P - pe - 2v 


3k 


U2\2 


3y 


(40) 


Dissipation rate c: 


- 

3x 


pv|£ 

3y 



y + 


^T ^ 
<?£ 3y 


r 3y 


+ C 






V p 


1/2 


( 41 ) 


with 


P - - 


-~i 1 dU — r I fdw w^ 

pu V 5^ - pw V - -) 


c F., P — 
y y e 


, - expf- 3 . 4 /(l + Rj,/ 50 )^] 


p - K 

T ve 


0 


k 


1 


a 


e 


1.3 


- 1.44 , C - 1.92 fl - .3 exp(-ld)] 

1 2 ^ 
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C "2 and pu'v 
^3 


„ ^ ag/r 

T 3y ’ ^ r " 3y 


(42) 


Dlacusslon about the coefficient will be presented later. In the preliminary 
calculatlona has been taken equal to 1. This coefficient represents the 
anisotropy existing between the Reynolds strers components and the mean strain 
components. ’n\ls model has been tested first, on an axlsyimmtrlc cylinder 
to check the effects of the corrections due to low Reynolds on the equation 
of dissipation rate e, through the coefficient We may see on Figure 1 

that the correction allows a 30 percent damping of the value of C. However 
this damping function starts to be effective at Reynolds numbers lower than 
two which means that occurs well Inside the viscous sublayer where the 
viscous effects are preponderant over the turbulent ones. In Figures 2 to 6 
we present some comparison between two calculations, the first one using the 
low Reynolds damping function In the equation of £ and the second one letting 
■ 1.92. It appears that the maximum differences around 1 percent on 
quantities such as the kinetic energy, the dissipation rate, the wall shear 
stress and below 1 percent for the velocity profj.le, which is well beyond 
the accuracy we may hope from this type of calculation and experimental 
measurements of these quantities. We may wonder why then lnq>lement such 
corrections. It Is possible In certain very low turbulence flows that this 
function may have some more Important effects. However it Is evident that 
again the changes due to the production terms In equation 40 and 41 are much 
more important and over a wider distance from the wall. Then It is In fact 
the modelling of the Reynolds stresses themselves which Is making the big 
difference . 


To illustrate this, the calculation of the boundary layer developing 
over a cylinder of which one part is rotating Is a good test case. In fact, 
the boundary layer first develops along the static cylinder as a classical 



26 


two dlnenslonal boundary layer. It is then aubmltted to a sudden transversal 
distortion by rotation of the downstream part of the cylinder. Sufficiently 
far away the origin of the spinning part, the reorganization of the three- 
dimensional boundary layer Into a two-dimensional one. In a relative coordinate 
system Is In an advanced phase. Very near the discontinuity In boundary 
conditions, the three dimensionality Is stronger. A recent paper of 
Arzoumanian et al. [25] shows that In a relative coordinate system linked to 
the moving cylinder, far from the origin of the spinning part, the turbulent 
quantities are very slightly Influenced by the three-dimensional effects, 
except the Reynolds stress uv. This suggests that the evolution of uv and wv 
in the fixed coordinate system are different and therefore the anisotropy of 
the Reynolds stress tensor Is quite Important, at least In the very beginning 
of the rotation. Then a model based on an Isotropic eddy viscosity (e.g. the 
two directions have the same characteristics) should fall to represent the 
dynamic characteristic of such a flow. 

We present hereafter some calculations of the rotating cylinder with an 
Isotropic eddy viscosity to Illustrate our purpose. The calculation Is 
confronted to the experimental results of Lohmann[201. Figures 7 and 8 show 
the developaient of the mean velocity profiles U and W, while Figure 9 shows 
the evolution of the limiting streamline angle along the streamwlse direction. 
Differences between the experimental and calculated profiles are up to 9 per- 
cent In our calculation. It can be easily seen In Figure 7. 

However, larger differences appear In the calculation of the Reynolds 
stresses u^^ and w'v^, (Figures 10 and 11) particularly for w*v* . It is then 
evident from these calculations that the Isotropic eddy viscosity falls to 
represent the characteristics of the flow. 


27 


Different authors have proposed sons alternative to oodel such a 
behavior. For example, Koosllln and Lockwood [21] derived eddy viscosities 
for the two directions and then define a viscosity ration 2 * 
which globally represent the anisotropy of the Reynolds stresses. Coustelx 
et al. [22] proposed a similar approach with the k-c model. 

A very simple way of deriving two different eddy viscosities is to 
use equation 38, which gives the ratio u^vVw*v* representing the differences 
between these two Reynolds stresses. 

Ue are now developing two different schemes using equation 38 to solve 
the turbulent field. The first scheme solves the momentum equations and the 
kinetic energy and dissipation rate equations first using the stresses u'v' 
andv^w^ calculated from the linear form of the equation 38 as we discussed 
previously. The second scheme solves the momentum equations and dissipation 
equation using again the stresses u*v* , v'w’ and u* , v* , w' calculated 
from the algebraic nonlinear system 38. The system of six nonlinear equations 
is solved with the iterative Newton-Raphson scheme. 

These two programs are now almost in their achieved form. A brief 
discussion of the numerical method used to solve the differential equations 


is presented in Appendix B. 



3. NUMERICAL AMALYSIS OF BLADE AND 
HUB WALL BOUNDARY LAYERS 


In the early stages of the present Investigation, an attempt was made 
to predict the blade boundary layer using the parabolized form of the 
Navier-Stokes equations. The procedure involves the solution of the momentum 
equations and the Poisson equation for the pressure. 

Numerous investigators have tried different ways of deriving and 
solving the Poisson equation for the pressure. In most of the cases the 
procedure requires some kind of iteration to ensure mass conservation. All 
the previous investigations indicated that the solution of the Poisson 
equation and the ensurance of mass conservation are the most troublesome 
parts of the whole procedure. The efficiency and accuracy of the method 
depends alisost exclusively on these factors. The situation becomes more 
difficult in the case of complex gecxmetry, like for example in the case of 
a turbomachinery rotor passage. 

Up ':o date, nobody has come up with a "clean" solution to the problem 
and all the evidence shows that coupling the momentum equation with the 
Poisson equation for the pressure is an ill-posed problem. Based on the 
above discussion we decided to try to solve the problem by using a procedure 
which does not involve the solution of a Poisson equation for the pressure. 

The first attempt was to solve the incompressible set of equations 
(three-TOmentum and continuity) with the vector of unknown given by. 


g • 


U 

V 

u 
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We tried to solve this system of equations by marching In the stream- 
vise direction C* We soon realized that: 

1. When treating the streamvlse pressure gradient fully explicitly 
(l.e., following the classical assumption of the parabolized Navler-Stokes 
equations) the Jacob Isn matrix, which expresses the changes In the streanwlse 
direction, becomes singular. 

2. When treating the streamvlse pressure gradient partly explicitly 
and partly Implicitly the solution Is unstable. This Is because our system 
of equations has both negative and positive eigenvalues. This makes It 
Impossible to solve the system of equations by marching In the streamwlse 
direction. 

An attempt Is presently being made to solve the compressible set of 
equations by marching In the streamwlse direction. This procedure was 
developed by Govlndan and Is described In ref. 24. The code Is presently 
working well, and It has already given very encouraging results for some 
test cases. So we decided to use that code for the solution of the flow In 
the blade passage. Including the hub wall and the blade boundary layer. 

We are presently working to adapt the code to our particular problem 
with the turbulence closure model described in section 2 of the report. 


4. EXPERIMENTAL STUDY OF THE ROTOR BLADE BOUNDARY 
LAYER IN AN AXIAL FLOW COMPRESSOR ROTOR 

The measurement of the blade boundary layer In the rotor of an axial 
flow compressor Is presently being carried out. The study is performed in 
the Axial Flow Compressor Facility, located in the Turbomachinery Labora- 
tory of the Department of Aerospace Engineering. The only simil4>r study 
available at the present time is that due to Lakshminarayana et al. [23,24] 
inside the rotor of an axial flow fan. 

All the measurements reported here are taken with a hot-wire probe 
rotating with the rotor. A miniature cross flow "X" hot-wire probe, 

TSI 1247, with sensor diameter of 3 urn and sensor length equal to 1 mt Is 
used. 

Both the sensors are within a circle of approximate diameter 1.5 mm. 

The sensors are located in the zR plane with their axis at 45* to the z axis 
(rig. 10). The probe is traversed normal to the blade surface. Since the 
flow traverse is done close to the blade surface, the component of velocity 
in the n direction Is assumed to be small. 

The boundary layer measurements are carried at non-dlmenslonallzed 
radii R - 0.583, 0.67, 0.75, 0,832, 0,918 and at various chordwlse locations. 
Only (me chordwlse location has been surveyed so far, at each of these 
radial locations. 

4 . 1 Experimental Results 

Only a brief Interpretation of the results is given here. Detailed 
Interpretation will be given when the whole set of data Is completed. 
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Th« itraaiBwlse velocity profiles ere shown In Figures 11 through 15. 

At all suction side locations the profiles are relatively well behaved. The 
profiles on the pressure side at the locations near to the hub (R ■ 0.583, 
0.67) ha>m an unconventional shape. This probably cones fron the secondary 
velocities developed In the hub wall region. 

The radial velocity profiles are shown In Figures 16 through 20. At 
the location nearest to the hub (R - 0.583) the radial velocity Is very snail 
(less than 4 percent of U^) . At R ■ 0.67 a strong radially outward velocity 
starts to develop very near the surface. At R > 0.75 we see that the radial 
velocity does not get stronger but It spreads away fron the surface (see 
Figures 18 and 19). At the location nearer to the tip, R - 0.918 (Figure 20), 
the radial velocity becones again small. Here the radial velocity Is affected 
by the leakage flow, and the annulus-wall botindary layer. 

The screamwlse and radial turbulent Intensities are shown In Figures 21 
through 26. In almost all locations we see an Increase of the turbulent 
Intensities when going towards the blade surface. 


5 . PUBLI CATIONS AND PRESENTATIONS 


The following papers were published during this period: 

1. Lakshmlnarayana, B., C. Hah, and T. R. Govlndan, "Three Dlnenslonal 

Turbulent Boundary Layer Development on a Fan Rotor Blade," AIAA 
Paper 82-1007, 1982. 

2. Pouagare, M., K. N. S. Murthy, and B. Lakshmlnarayana, "Three Dimensional 

Flow Field Inside the Passage of a Low Speed Axial Flow Compressor 
Rotor," AIAA Paper 82-1006, 1982. 

The following presentations were made during the reporting period. 

1. B. Lakshmlnarayana, "Experimental Study of the Rotor Blade Boundary 

Layer in a Single Stage Fan," lUTAM Syiig)oaium on Three Dimensional 
Turbulent Boundary Layers, Technlsche Unlversitat Berlin, Germany, 

March 29, 1982. 

2. B. Lakshmlnarayana, "Three Dimensional Turbul.int Boundary Layer Develop- 

ment on a Fan Rotor Blade," AIAA/ASME Third Joint Thermophysics, 

Fluids, Plasma and Heat Transfer Conference, St. Louis, Missouri, 

June 11, 1982. 

3. M. Pouagare, "Three Dimensional Flow Field Inside the Passage of a Low 

Speed Axial Flow Compressor Rotor," AIAA/ASME Third Joint Thermophysics, 
Fluids, Plasma and Heat Transfer Conference, St. Louis, Missouri, 

June 11, 1982. 
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D«rivaclon of 
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b«caus* the derivation at point x and are independent. It la atvaight- 
forvard in the case of carteaian coordlnatea. In the caae of generalised 
tenaor see Arls (12]. let ua define new independent variablea aa follows: 
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Then with the hypothesis of hosngcneous flow the U,^„ is independent of the 
point (I), therefore It can be shown [3,10,14] 
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m the Indices n* and r' are dusny se c«i rewrite equation A*2 as follows: 
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Patankar-Spalding Numerical Method 
for Two-Dimensional Botmdary Layers 


1. Intrjd action 

We are gJing to give some indications on the numerical method established 
by Patankar and Spalding in 1967 and modified in 1970 by these authors. 

The presentation will be brief and for further information the reader may 
want to refer tr r e b:>ok entitled HEAT AND MASS TRANSFER IN BOUNDARY LAYERS— 
A GENERAL CALCULATION PROCEDURE, by S. V. Patankar and D. B. Spalding, 
Intertextbook, London, 1970. The method is of the marching type and is 
adapted for two-dimensional boundary layers. The boundary layers approxima- 
tions lead to a system of parabolic, partial derivatives equations. With the 
Patankar-Spalding method we may calculate a purely dynamic boundary layer, 
as well as a thermal boundary layer, or whatever contaminant evolution. 

This method does not need any boundary layer similitude hypothesis, 
and the principal characteristics of the method are the following: 

— A non-dimensional streamfunction is used as independent variable, 
then the couple (x,y) of the physical plane is replaced by the couple (x,w) 
which allows an automatic adaptation of the thickness at each abscissa x. 

— The implicit scheme of discretization using the micro- integral 
technique, verifies the conservative properties of the partial derivatives 
equations and is unconditionally stable. 

— The marching solution is very well adapted to calculate a flow, 
very quickly and using a really few memory storage in the computer. 

— To invert the matrices, the method uses a tridiagonal algorithm 
derived from the Gauss method, this is done witn \ . simple iterative 


formulas . 
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2. Formulation of the Equation 

Here, ve will use a tvo-dimenslonal cartesian frame to present the 
derivation of the equations. 

For each transport quantity ^ of Interest In the botindary layer. Its 
equatl<« can take the form: 


P U 


3<j) 

3x 


+ p V 


30 

3y 


3r 


eff 


30/ 3y 


3y 


+ S. 


( 1 ) 


where Is an effective diffusivlty coefficient and represents a 

"source” term. 

The continuity equation is joined to equation (1) . 

Using Von Mises transformation, where (x,y) is replaced by (x,0) , 0 
is the streamf unction (p U * 30/3y), equation (1) becomes; 


30 ^ ^eff ^0 

3x ” 30 PU 


( 2 ) 


In this new coordinate system (x,0), the grid is automatically adapting to 
the growth of the boundary layer thickness. However, for simplicity in the 
definition of boundary conditions, it is useful to introduce another coordinate 
system (x,w) where w is defined by w » 0 - 

respectively the values of 0 at each boundary. We have then; 


w. 


w * 1 
E 


d0. 


d0. 


dx ’ dx 

In this new coordinate system, equation (1) becomes:’ 


-P V 
^E E 


30 , . 30 3c 30/ 3w , 

+ (a + bw) » ;; + d 


3x 


3w 


3w 
dif , 


0 

source 


( 3 ) 


conv. 
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where a and b are functions of x 


♦e'*i 


°E ^'e ' ‘*1 

♦j- »I 


where c and d. are functions of x, w, ♦ 
9 


pur 


eff 


(♦e ■ * 1 ^' 


d^ depends on the definition of the quantity 9 
(for example for velocity U, d . is the pressure 


gradient) 


9 


The mass transfer rates m" * pV through the boundaries 1 and E are specified 

by the nature of the boundaries (symmetry axis, wall, free boundary, •••). 

In the case of a symmetry axis or a wall, these rates are zero. Eor a free 

boundary G, Patankar and Spalding express the mass transfer rate as a function 

of the effective diffusion, as follows: 

3/3y (y^jf 3U/3y)' 

3u73y 

This necessitates the knowledge of the effective viscosity. We may note that 
it is this transfer rate which controls the "entrainment" of the marching 
solution, then it is important to control the definition of this mass transfer 
rate m". Some damping functions may be used to diminish the exchange of flow 
through the boundary G when that is necessary (see Patankar-Spalding) . 

Now, equation (3) is a non-linear partial derivatives equation with 
boundary conditions depending on the specificity of the physical problem. Then, 
to solve this equation, in the general case, we have to use a finite difference 
technique which is described hereafter. 


m^ “ lim 
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3. Numerical Discretization 

Equation (3) la discretized and Integrated by a marching method In 
direction x. At each step of the Integration* the values of ^ are known at 
each point of the grid in direction w for a given section* and we may get 
the values of ^ for the next section. 

Generally, to obtain a finite difference form of a partial derivative 
equation we may use the Taylor series. But It Is also possible to obtain 
the finite difference equation considering each term of the non-dlscretlzed 
equation as an average Integration on a smaj 1 • ume around each point of 
discretization; this volume Is called the control voltane. This Is the 
process which is used by Patankar-Spaldlng, Introducing certain hypotheses 
relatively to the nature and the variation of 0 between grid nods. Two 
assunqjtlons are made on the variation of First It Is supposed* that 4> 
varies linearly between nodes w^, and secondly, 0 is set constant In the 
interval Ix^, x^l and equal to Its value In x^. Then, for the non-boundary 
zone, equation (3) may be Integrated In the rectangular dot domain presented 
below (Figure 1) as follows: 


+ + 

[/ (})j^dw - / (Jiydw]/(rx+ [{(a + + bw)<{)jj}_] 


+ 

- b / 4»jjdw 


'(c ii] . fc ' 


+ 

+ / d ,dw 


( 4 ) 


Equation (4) can be put, then. In the recurrent form: 


(t ^ A 0 +B4> +C 

^Di i ^ °i-i ^ 


( 5 ) 


where A^, are coefficients and are known at each node. 


40 



Figure 1. 


In the case of the boundary region, a different scheme is used. On 
boundary I, for example (Fig. 2) a fictitious value called the "slip 
value" is introduced such as the linear interpolation of ()> between nodes W 2 
and gives the right slope at the intermediate point «2 In fact, the 
boundary regions are important as they control the fluxes and entrainment 
rate, the variations of 4 >(w) could not be supposed linear in a domain which 
is only half of an equivslent domain away from the boundaries. The values 
of 4> at the boundaries 1 and E are given by the boundary conditions; if 
the grid is numbered from 1 to N + 3, we need two discretized equations to 
relate and to, respectively, and 0^, and and Then 

we may write; 

^2 “ ^ 2'^3 ® 2*^1 ^2 ’ ‘** N +2 “ N + 2 '^ N +3 ® N + 2 '^ N +1 S +2 

where A 2 , 82 * C 2 and ’ ®N+2* ^N+2 known and include the boundary 


conditions . 
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Figure 2. 


Then, we obtain the following tridiagonal system: 


"**1 “ Vi+1 ®A-1 ^i - 2, .... N+2) 


( 6 ) 


The system is solved, using the so-called Thomas algorithm: 


“ ®i‘*’i+l ^i - 2. .... N+2) 


with 


^ ■ ^i^^i-l 


°i ^i-1 ^ ^i 
^ ~ ^i^^i-l 


(i - 2, .... N+1) (i - 2 N+2) 


^1 " ^1 




Then the values 4>(x,w) are known for each w; if we want to locate the value 
in the physical plane (x,y), we have to compute the y^ at each station x. 

In the particular case of a plane flow, the value y corresponding to w is 


^i 



J PU 


given by; 
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Entrainment Rate 

In the case of a wall at the boundary 1(E), la given by the 

wall conditions (blowing, suction, etc.)* 

In the case of a free boundary, Patankar and Spalding have proposed the 
following formulations: 



^eff 


2.5 


73 - 7j 


eff, 


N+1.5 


^N+2 “ ^N+1 


which come from the expression: 


.. ,, . , 3^0, 8b] 

5 JZ ~ 

_ 9y j 


where the second term Is neglected. Then the entrainment rate is always 
positive, and by the way, instabilities of the kind positive-negative- 
positive in calculation of the rate are avoided. 


Source l€ ms 


Different procedures are applicable to discretize the sources S 




i+1/2 

S* ■ ('J'r " 

S . is linearized in <p and is written as; 




For example, for the U component equation the source tern is the 


pressure gradient, and is given by: 
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i 




^ 1+1 


2 





0 


In the general case, for exan^le, total enthalpy, kinetic energy, 
dissipation rate, many decompositions are possible. 

We may put the positive sources in and the negative ones in S^, or 
put the proportional sources to i° while the others are put in S^, 
or also put all the terms in and in divided by (Launder) . In fact, 
it does not seem that a particular technique is better than another one. 
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CH^^RT 












XM0.64 CM 

Figure 3. Coaparlson of calculation with and without function 
Turbulent kinetic energy, k/U^ 








Figure 6. Conqparlson of calculation with and without function 
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Figure 7. Comparison with Lehman 's data; W^/Ug - 1.45 (where la the speed of rotation of the cylinder) 
Axial velocity profile 





53 


ORIGIP^AI. 13 



J33*. 





Figure 11. Comparison with Lohman’s data; Wc/Ug “ 1.45 (where Wg is the speed of rotation of the cylinder) 
Shear stress In tangential direction 
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Figure 12. Nature of blade boundary layer and 
notations used 
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Figure 15. Streanwlse velocity (U^) profiles, R - 0.75 
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v’ a: uj < z: 3 H (/) uj >uj_jouhh-> 



Figure 16. Streaiwlse velocity (U^) profile, R - 0.75, * - 0.97 (P.s. 
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</) I- a uj < z :» H <0 uj >uj_iouhi-> 


Figure 17. Streamlse velocity (U ) profile, R ” 0.918 
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Figure 18. Radial velocity (W) profiles, R - 0.583, e - 0.68 (P.S.), z - 0.84 (S.S.) 
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Figure 19. Radial velocity (W) profiles, R - 0,67, z - 0.68 (P.S.), z - 0.86 (S.S.) 
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Figure 23. Streamwise turbulent intensity (x ) profiles, R » 0.75, z = 0.69 (P.S.),z - 0.88 (S. 
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Figure 25. Streamwlse turbulent intensity (t ) profile, R * 0.918, z * 0.68 (P. 
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Figure 26. Radial turbulent Intensity profiles R “ 0.75, z * 0.69 (P.S.), z » 0.88 (S.S. 
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Figure 27. Radial turbulent intensity profile (t_)» R “ 0.75, z - 0.97 (P.S.) 
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Figure 28. Radial turbulent intensity profiles (t^), R “ 0.918, z - 0.68 (P.S 


